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ABSTRACT 

' Using a sample of 2408 time-resolved spectra for 91 BATSE GRBs presented by Preece et al., we 

, show that the relation between the isotropic-equivalent luminosity (iiso) and the peak energy {E^) 

^ \ of the vFi, spectrum in the cosmological rest frame, Ljso oc E^, holds within these bursts, and also 

■ holds among these GRBs, assuming that the burst rate as a function of redshift is proportional to the 
star formation rate. The possible implications of this relation for the fireball models are discussed by 
defining a parameter lo = (Liso/lO^^erg s~^)°'^/(i?p/200 keV). It is found that lo is narrowly clustered in 

' 0.1 — 1. We constrain some parameters for both the internal shock and external shock models from the 

requirement of a; 0.1 — 1, assuming that these model parameters are uncorrelated. The distributions 
of the parameters suggest that if the prompt gamma-rays are produced from kinetic-energy-dominated 
internal shocks, they may be radiated from a region around R ^ 10^^ — 10^^ cm (or Lorentz factor 
~ 130 — 410) with a combined internal shock parameter ^ 0.1 — 1 during the prompt gamma-ray 
phase, which are consistent with the standard internal shock model; if the prompt gamma-rays of 
these GRBs are radiated from magnetic-dissipation-dominated external shocks, the narrow cluster of uj 
^ . requires cr - 1 - 470, F - 216 - 511, E - 10^^ - lO^"* ergs, n - 0.5 - 470 cm'^, and Ce ~ 0.36 - 3.6, 

' where a is the ratio of the cold-to-hot luminosity components, F the bulk Lorentz factor of the fireball, 

, E the total energy release in gamma-ray band, n the medium number density, and a combined 

■ external shock parameter, which are also in a good agreement with the fittings to the afterglow data. 
' These results indicate that both the kinetic-energy-dominated internal shock model and the magnetic- 

pi . dissipation-dominated external shock model can well interpret the L\so oc E^ relation and the value of 
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. 1. INTRODUCTION population of X-ray hashes, which are thought to be a 



n u 4- rr^-DTD \ u 1- J i u low energy extension of typical GRBs (Heise et al. 2001, 

Gamma-ray bursts (GRBs) are now believed to be pro- , . „„„„n i .i • i ;• a.. • 

J J , . , J u J- 1 ■ -iu 1 Kippen et al. 2003). Based on this relation, Atteia 2003) 

duced by lets powered by central engines with a stan- , ^ ^ , , ' . , , , -n • i- . r 



1-1 j-j- { also constructed a simple redshift indicator for GRBs. 

. , dard energy reservoir at cosmological distances (see se- „ , . \. , ■ , ;• i i i ■ 

?H • • u T?- u p lAAc n- innn One may ask: whether or not this relation holds m 

■ ries reviews by irishman & Meegan 1995; Piran 1999; van ... ^t^t^o mi 

■ ■ ' Paradijs et al. 2000; Cheng & Lu 2001; Meszaros 2002; f^^^ segment withm a GRB? The answer remains un- 

Zhang & Meszaros 2003). ^^^^^"^ P^^^*^^^' combining the results 

rri, i • • r i c r^TDo j-u j- j- mentioned above, one might suggest that this relation is 

1 he most impressive features of GRBs are the great di- . i , i • ; i ? , ^ 

r ,1 . ,. , , J i 1 1. 1 • J a universal law during the prompt gamma-ray phase, and 

versities of their hght curves and spectral behaviors, and , , ■ , r. i ii i i t .i ■ t , 

, 1 , , . •, ■ , 11 n, presents some constraints on fireball models. In this Let- 
extremely large luminosities. Ihese spectra are well fit- . ,. ,i . . tt • ^ mA^n 
ted by the Band function (Band et al. 1993). However, t^"^' we investigate this issue. Using a sample of 2408 time- 
.1 1- ,• 1 • J- 1 J ■ j-i, +1, resolved spectra ior 91 BAiSE GRBs presented by Preece 
the radiation mechanism at work during the prompt phase , , ,„„fls , , , , ii .,, . ,i 

1 J 4- J Au^, ^, i- ^ u et al. (2000), we show that this relation holds withm these 

remains poorly understood. Although the spectral be- , , ^ /', , , , ,, ^T^T^ • .i , 

1 . 1 1 ■ -J. J n j-jx i- c bursts, and also holds among these GRBs, assuming that 

havior and the luminosity are dramatically ditterent from , , ' , , ^ , . ^ ^ i i • n. • . • i , 

1 J. J.-U ■ -L ■ -lii • -4. T the burst rate as a lunction ot redshiit is proportional to 

burst to burst, the isotropic-equivalent luminosity, Lisa, i , r , ttt , i i i i i • 

/ . -14. j-4.jT,4.u the star formation rate. We suggest that both the kmetic- 

(or isotropic-equivalent energy radiated by the source, , . ; i . ; i i r i i i i 

. / • r n -1 energy-dominated internal shock model and the magnetic- 

i^i^o), and Ep, the peak energies of i^F, spectrum m the dissipation-dominated external shock model may well in- 

rest frame among GRBs, obey an empirical relation ot terpret this relation 

Lisa oc Ep^ (Amati et al. 2002; Yonetoku et al. 2003; Throughout this work we adopt iJo = 65 km s^^ Mpc-\ 

Sakamoto et al. 2004; Lamb et al. 2003a, b, c). This fi^^ = 0.3, and f^A = 0.7. 
relation was revisited in standard synchrotron/inverse- 

Compton/synchro-Compton models (Zhang & Meszaros 2. THE Liso - Ep RELATION WITHIN A GRB 

2002). Recently, Sakamoto et al. (2004) and Lamb et i ^ i i i , ' 

al. (2003a, b, c) pointed out that HETE-2 observations ^ithm a GRB, the relation between and E^ is 

not only confirm this correlation, but also extend it to the equivalent to a relation between the observed flux {F) and 

peak energy, Ep. Thus, we examine whether or not both 
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F and Ep follow a relation of F (x E^. The time-resolved 
spectral catalog presented by Preece et al. (2000) includes 
156 long, bright GRBs. Four spectral models were used 
in their spectral fittings. Different models might present 
different fitting results. Among 156 GRBs, 91 GRBs were 
fitted by the Band function (Band et al. 1993). We only 
include these GRBs into our analysis. There are 2408 
time-resolved spectra for these GRBs. In our analysis, the 
values of Ep are taken from this catalog. The data used 
for spectral fittings were observed by different BATSE de- 
tectors. Nominal energy coverage is 25 to 1800 keV, with 
some small variations between these detectors. The fluxes 
presented in the catalog are in an energy band correspond- 
ing to the detectors. Hence, the values of F in our analysis 
are not the fluxes presented in the catalog, but are the inte- 
grated fluxes in energy band 30-10000 keV (e.g., Yonetoku 
et al. 2003) derived from the model spectral parameters 
in the catalog. 

We evaluate the relation of F (x E^ within a GRB by 
the linear correlation coefficient of the two quantities, log 
F and log E^. We calculate the linear correlation coefH- 
cients (r) and chance probabilities (p) for each burst with 
the Spearman rank correlation analysis method. The dis- 
tribution of r is presented in Figure 1. We flnd that about 
75% GRBs exhibit a strong correlation between the two 
quantities with r > 0.5 and p < 0.0001. We illustrate 16 
cases in Figure 2. These results show that this relation 
holds within these GRBs. 

To examine whether or not this relation holds among 
these GRBs, we assume that the redshift distribution for 
these GRBs is same as that presented by Bloom (2003). 
Bloom (2003) assumed the burst rate as a function of red- 
shift is proportional to the star formation rate as a function 
of redshift, and presented the observed redshift distribu- 
tion incorporating with observational biases (SFRl model 
is used in this work, see Porciani & Madau 2001). We 
derive a value of redshift for a given GRB from this distri- 
bution by a simple Monte Carlo simulation method. To do 
so, we flrst derive the accumulated probability distribution 
of the Bloom's redshift distribution, P{z) (0 < P{z) < 1), 
then generate a random number for a given GRB, m 
(0 < m < 1), and flnally obtain the value of z from 
P{z) = m, i.e., z = P~^{m). Hence, we calculate the 
values of -Liso,52 and Ep{l + z) for these GRBs, where 
-^^iso,52 is in units of 10^^ erg s^^. The iiso,52 as a function 
of £'p(l + z) is shown in Figure 3. The linear correla- 
tion coefficient of the two quantities is 0.63 with a chance 
probability p < 0.0001. A best fit with a model of 
^iso,52 [ii'p(l + z)]^ is shown in Figure 3 (solid line). 
The reduced is 0.23. These results well suggest that 
the relation of iiso,52 oc Ep holds among these GRBs. 

3. IMPLICATIONS FOR FIREBALL MODELS 

The above results well suggest that the relation of Liso oc 
Ep remains within a GRB. This implies that the relation 
is independent of the temporal evolution of a fireball. This 



e~^, see Zhang & Meszaros (2002). 
2Wc check Eq. (23) in Zhang & Meszaros (2002) 
In order to discuss conveniently the distribution of a, 
E^Qiigh a) ~ 2m\<ieYCe'T\'^Tl'^L\{^E-^'^n\'^{l + z)- 



might provide strong constraints on fireball models. We 
define a quantity, uo, to discuss these possible constraints, 
which is given by 



LO = 



rl/2 

^iso,52 



E, 



(1) 



'p,200 



where Ep is the peak energy in the rest frame in units of 
200 keV. Since Lobs = S^L^o and E^hs = S^cot we obtain 
'^obs = Wco, where S is the Doppler-boosting factor, and 
the subscripts "obs" and "co" denote the observer frame 
and the comoving frame, respectively. This indicates that 
Lo is not influenced by the Doppler-boosting effect, imply- 
ing that it seems to be an intrinsic parameter relevant to 
the properties of the fireball models. 

Two competing fireball models involve the internal 
shock and external shock models. Zhang & Meszaros 
(2002) analyzed correlations between Ep and other pa- 
rameters in several cases. From their results, we find 
that the kinetic-energy-dominated internal shock model 
(low cj, where a is the ratio of the cold-to-hot luminosity 
components) and magnetic-dissipation-dominated exter- 
nal shock model (high a) have a potential to interpret the 
I/iso (X E^ relation. Prom their Eqs. (17) and (18), we 
derive 

w ^ CiRi3 (2) 

for the internal shock model with low a, where d is a com- 
bined internal shock parameter^, and R13 the radius of the 
fireball in units of 10^^ cm. Prom their Eq (23)^, we obtain 



l/3y^-8/3 pl/3 -1/3 



o.36Cc'Tr'''r2-°'^S5,^ 



(3) 



for the external shock model with high a. 

From Eq. (2) one can see that, for the internal shock 
model, CO is determined by the radius of the gamma- 
ray emission region (hence the bulk Lorentz factor since 
i? oc F^) and ([. Please note that d is related to the 
shock parameters, such as the index of the electron dis- 
tribution, the electron equipartition factor, the magnetic 
equipartition factor, the pitch angle of electrons, the rel- 
ative Lorentz factor between the shells, etc. These pa- 
rameters arc only related to the physic;s of colliding shells. 
However, for the external shock model, the case becomes 
more complicated. Eq. (3) shows that ui is determined by 
the parameters of both the shocks and surroundings. 

We calculate uj for each temporal segment within a GRB 
and investigate the evolution of w by a linear model of 
u) oc kt. The value of k evaluates the general trend of 
the temporal evolution feature of oj: the larger the abso- 
lute value of k is, the more significantly uj evolves. The 
distributions of lo and k are shown in Figures 4 and 5, re- 
spectively. We find that oo mainly distributes in 0.1 — 1, 
and k is narrowly clustered in -0.03 — 0.03. These results 
show that to seems to be an invariant without temporal 
evolution for different GRBs and even for different tempo- 
ral segments within a GRB . 



and find that the coefficient of this equation is 12 keV, rather than 880 keV. 
we also assume a^^^/{l + a)^^^ ~ cr^/"^. Thus, the equation is re-scaled as 
where Ce is a combined external shock parameter (^e = ^^^-1^7 ^ sin~^ ^f, see 



Zhang & Meszaros 2002), ai 
of 10 cm-3. 



■ O-/10, r2.5 = r/lO^-S, E53 the total energy in units of 10 ergs, and ni the medium number density in units 
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The parameters in Eqs. (2) and (3) seem to be un- 
correlated, although we do not know if this is really 
the case. We simply assume that they are uncorrelated, 
and constrain their distributions from the requirement of 
a; 0.1 — 1. We suggest that these parameters should be 
clustered in the same range as that of lo. Thus, we de- 
rive R - 10^2 - lO^-' c;m and Ci 0.1 - 1 for the internal 
shock model, implying that most of the gamma-rays are 
radiated from a region around R ~ 10^^ — 10^^ cm with 
similar shock-acceleration and radiation mechanisms dur- 
ing the prompt gamma- ray phase. Since R ~ 2r^C(5iv ~ 
0.6 X 10-^^r| 5(5iv,-3cm , where c is the speed of light, and 
fv,-3 the variability timescale in units of 10~^ second, we 
obtain F ~ 130 — 410. These parameters are consistent 
with the standard internal shock model. For the exter- 
nal shock model, we derive a ~ 1 — 470, F ~ 216 — 511, 
Ce ~ 0.36 - 3.6, E ~ 10^1 - 10^^ ergs, and n ~ 0.5 - 470 
cm~^. The distributions of these parameters are in a good 
agreement with the fittings to the afterglows (Panaitescu 
& Kumar 2001). 

The above results indicate that both the low a internal 
shock model and the high a external shock model can well 
interpret the L oc relation and the value of w. 

4. CONCLUSIONS AND DISCUSSION 

Using a sample of 2408 time-resolved spectra for 91 

long, bright GRBs presented by Preecc et al. (2000), 
we show that the L\^o ^ relation holds within these 
BATSE bursts, and this relation also holds among these 
GRBs by assuming that the burst rate as a function of 
redshift is proportional to the star formation rate. 

We discuss possible implications of this relationship 
for the fireball models by defining a parameter a; = 



(iiso/10^2^rg s-i)0-V(£^p/200kcV). It is found that w is 
not influenced by the Doppler boosting effect, and it is de- 
termined by the gamma-ray emission region and shock pa- 
rameters in the kinetic-energy-dominated internal shock 
model or determined by the parameters of both the 
shock and the environment in the magnetic-dissipation- 
dominated external shock model. We derive the distri- 
butions of some parameters for both the internal shock 
model and the external shock model from the requirement 
of a; ^ 0.1 — 1. We suggest that if the prompt gamma- 
rays are produced from a kinetic-energy-dominated in- 
ternal shock, they may be radiated from a region around 
R ~ 10^2 _ iQi3 pjjj (qj. Lorentz factor ~ 130 - 410) with 
an internal shock parameter Ci ~ 0.1 — 1, which is consis- 
tent with the standard internal shock model; if the prompt 
gamma-rays of these GRBs are radiated from magnetic- 
dissipation-dominated external shocks, the a; ~ 0.1 — 1 
requires cr - 1 - 470, F ~ 216 - 511, Co 0.36 - 3.6, 
E - lO-^i - 10''^'' ergs, and n - 0.5-470 cm--\ Please note 
that the distributions for these model parameters for both 
the internal and external shock models are based on the 
assumption that they are uncorrelated. Although these 
parameters seem to be uncorrelated, we do not know if it 
is really the case. If these parameters are correlated dur- 
ing prompt gamma-ray phase, these distributions are not 
valid. 
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Fig. 1. — Distribution of the linear coefficients for log_F-logi7p . 
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Fig. 2. — The observed flux as a function of Ep for 16 GRBs. 
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Fig. 3.— The Liso,B2 as a function of Sp(i+^) for 2408 GRB spectra. The sohd Une is Liso,52 = IQ-^ i x [Ep{l + z)]^. 



7 




Fig. 4. — Distribution of w for 2408 spectra. 



8 




0.25 



Fig. 5.— Distribution of k for 91 GRBs. 



